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ABSTRACT 

The atmospheric parameters of the components of the 16 Cygni binary system, in which the sec- 
ondary has a gas giant planet detected, are measured accurately using high quahty observational data. 
Abundances relative to solar are obtained for 25 elements with a mean error of (t([X/II]) = 0.023 dex. 
The fact that 16CygA has about four times more lithium than 16CygB is normal considering the 
slightly different masses of the stars. The abundance patterns of 16CygA and B, relative to iron, 
are typical of that observed in most of the so-called solar twin stars, with the exception of the heavy 
elements {Z > 30), which can, however, be explained by Galactic chemical evolution. Differential 
(A-B) abundances are measured with even higher precision (cr(A[X/II]) — O.OlSdex, on average). 
We find that 16CygA is more metal-rich than 16CygB by A[M/H] = +0.041 ± 0.007 dex. On an 
element-to-element basis, no correlation between the A-B abundance differences and dust condensa- 
tion temperature (Tq) is detected. Based on these results, we conclude that if the process of planet 
formation around 16CygB is responsible for the observed abundance pattern, the formation of gas 
giants produces a constant downwards shift in the photospheric abundance of metals, without a Tq 
correlation. The latter would be produced by the formation of terrestrial planets instead, as suggested 
by other recent works on precise elemental abundances. Nevertheless, a scenario consistent with these 
observations requires the convective envelopes of ~ ^Mq stars to reach their present-day sizes about 
three times quicker than predicted by standard stellar evolution models. 

Subject headings: stars: abundances — stars: fundamental parameters — stars: planetary systems — 
stars: individual (16CygA, 16CygB) 



1. INTRODUCTION 

It is often assumed that the members of multiple stel- 
lar systems with similar components have the exact same 
chemical composition since they formed from a common 
molecular cloud and have undergone similar evolutionary 
paths. A number of studies have found, however, that 
elemental abund ance differences in binar y stars are not 
unco r rimon (e.g.. iGratton et al.l: 2001: La ws fc Gqnzalea 
[200lt ISadakane et all 120031: IDesidcra ct al.l 12004 I2'006ir 
Systematic errors in the spectroscopic analysis could be 
blamed for some of these anomalies but small differences 
have been observed even in pairs with very similar com- 
ponents, i.e., in cases where differential analysis is ex- 
pected to minimize the errors. The origin of these abun- 
dance anomalies remains unknown. 

Stars and their planets form nearly at the same time. 
Therefore, it is likely that the process of planet forma- 
tion leaves detectable signatures in the chemical composi- 
tion of the host stars. The well-known planet-metallicity 
correlation in late-type main-sequence stars, for exam- 
ple, suggests that high metallicity increases the proba- 
bility of forming planets, at least the gas giants (e.g., 
Marcv et a l.' '2005: 'Udrv fc Santos' '20071 : iJohnson et alj 
2010D iMelcndcz ct al. (200_a, hereafter M09) have fur- 
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ther suggested that small anomalies detected in the so- 
lar chemical composition compared to that of the ma- 
jority of so-called solar twin stars, objects with spectra 
nearly indistinguishable from the solar spectrum, could 
be due to the forni ation of terrestrial planets (see also 
iRamirez et al.|[2"009l hereafter R09). They detect a small 
deficiency of refractory elements in the solar photosphere 
which would disappear if the rocky material that once 
formed in the solar system were diluted into the present- 
day solar convective envelope (|Chambersll2010( ). Thus, 
the missing mass of refractories in the Sun seems to have 
been locked-up in the rocky bodies of the solar system, 
including the terrestrial planets. 

The hypothesis described above assumes that the for- 
mation of the massive gas giants in the solar system did 
not affect the solar chemical composition, which seems 
unlikely given the much larger masses of these planets 
in comparison with the terrestrial planets. Nevertheless, 
this could be the case if, for example, the abundance ra- 
tios of different metals in the gas giants combined were 
similar to those found in the Sun. Unfortunately, the 
detailed chemical composition of the gas giants, in par- 
ticular their cores, is not yet known with enough pre- 
cision to rule out or confirm this assumption. Careful 
spectroscopic studies of stars with gas giant planets de- 
tected could therefore provide us with important clues 
about the impact of planet formation on the chemical 
composition of the host stars. 

The binary star 16 Cygni is an ideal system to look 
for small elemental abundance anomalies and to study 
the impact of gas giant planet formation. 16 Cyg has 
long been known as the nearest bound pair of solar twin 
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stars (e.g.. iFriel et all 119931: JPorto de Mello fc da Silval 
Il997t iHauser fc Marcvlll999[ ). although strictly speak- 
ing they should be referr ed to as solar analogs (cf. 
iMelendez fc Ramired 120071 ) . in particular because they 
appear to be older than the Sun. Careful line-by-line 
differential analysis of this type of stars yields elemen- 
tal abundances with the highest precision possible, as 
demonstrated by the works of M09 and R09. More- 
over, the components of the 16 Cyg system are bright 
(y ~ 6 mag) and therefore very high quality spectra can 
be obtained for them using mid-sized telescopes. Thus, 
extremely high precision differential spectroscopic anal- 
ysis can be applied to this system. 

The secondary star in the 16 Cyg binary system, 
16 Cyg B, is known to host a gas giant planet with a 
minimum mass of 1.5 Jup iter masses (1.5 Mj) i n a highly 
eccentric (e ~ 0.6) orbit. iCochran et al.l (|1997[ ) were the 
first to detect this planet and they also obtained accu- 
rate radial velocity data that showed that 16 Cyg A does 
not host the same type of planet (gas giant in short- 
period orbit). Up-to-date follow-up observations of this 
pair of stars have not changed their planetary status 
(W. Cochran, private communication). The presence of 
low-mass and/or long orbital period planets around ci- 
ther of these two stars is, of course, not to be ruled out. 

A faint companion, likely a cool M-dwarf, to the pri- 
mary star 16 Cyg A, at an angular distance of about 
3.4 arcsec has been reported by a number of authors (e.g. 



Hauser fc M arcv°1999^' Turner et al.ll2001l : lPatience et alJ 
2002). The absolute magnitude difference between 



16 Cyg A and C (the faint companion) is about 8.6 mag in 
the R band. Their reported separation is large enough 
for the light from 16 Cyg C to not interfere with that 
from 16 Cyg A in our data but even if 16 Cyg C was su- 
perposed with 16 Cyg A at the time of our observation, it 
would affect the combined spectrum by less than 0.04 % 
at A = 6000 A and even less at shorter wavelengths. 

The atmospheric parameters of the 16 Cyg pair have 
been measured by a number of authors using a variety 
of techniques. Few of these previous studies, however, 
present results based on extremely high precision anal- 
ysis, i.e., significantly better t han "standard" rnethods . 
One of these studies is that bv lLaws fc GonzaIe3 ()200H ). 
who find that 16 Cyg A is more iron-rich than 16 Cyg B 
by about 6 ± 2%. They also comment as a note added 
in proof that a significant positive correlation with dust 
condensation temperature is found for the abundance dif- 
ference of 13 elements analyzed. They claim that their re- 
sults support a self-pollution model in which 16 Cyg A ac- 
creted planetary material that once formed there. They 
also argue that this scenario could explain the observed 
differ ence in lithium ab undance between 16 Cyg A and B 
(e.g.. iKing et al.lll997D . Later work bv iTakedal (|2005[ ). 
however, suggests that 16 Cyg A and B have exactly the 
same chemical composit ion and attributed the discrep- 
ancy with the results by iLaws fc Gonzalea to differences 
in the adopted stellar parameters. 

Here we derive atmospheric parameters and elemental 
abundances as precisely and accurately as possible using 
high quality data of 16 Cyg A and B. High precision ele- 
mental abundances are presented for 25 elements in the 
16 Cyg binary system. We discuss scenarios that could 
explain the observed abundance patterns. 



2. DATA 

The photometric data used in this paper, i.e., mag- 
nitudes and colors, correspond to those listed in 
the General Catalogu e of Photometric Data (GCPD, 
iMermilliod et al.l 1 199 71) . We adopted the mean values 
(and errors) given in the GCPD, which have been com- 
puted using data from several different sources. For the 
(V, B — V) pair, for example, more than 40 observations 
from about 10 sources are available for both 16 Cyg A 
and B. 

To calculate the distance to the 16 Cyg system, we 
adopted the measured trigonometri c parallaxes from th e 
latest reduction of Hipparcos data (|van LeeuwenI 120071 ). 
Averaging the published values for the two stars we de- 
rive a parallax of 47.29 ± 0.21 mas, which corresponds to 
a distance of 21.15 ± 0.09 pc. 

High resolution {R = A/AA ~ 60,000), high signal-to- 
noise ratio {S/N ~ 300 — 500, per pixel) spectra were ob- 
tained by us with the R. G. Tull coude spectrograph on 
the 2.7m Harlan J.Smith Telescope at McDonald Ob- 
servatory on April 25, 2011. These spectra cover the 
wavelength region from 4130 to 10 000 A, complete up 
to about 5935 A, and with increasingly wider gaps of 2 
to 120 A towards redder wavelengths. The spectral win- 
dows from 6360 to 6440 A and from 6480 to 6520 A are af- 
fected by the so-called "picket fence," a straylight source 
of e mission lines wit hin the spectrograph (see Sect. 4.6 
in iTull et all 119951 ) , and therefore were not used in 
this work. Two exposures of 20 minutes each were ac- 
quired for both 16 Cyg A and B, one immediately after 
the other. Data were reduced in the standard manner 
with the IRAF echelle package, using identical reduc- 
tion parameters for the two stars (e.g., polynomial orders 
for the overscan and scattered light removal as well as or- 
der tracing and aperture widths, fiat fields, wavelength 
solution, etc.). The spectrograph was set so that the Ha 
line at 6562 A fell in the center of the detector. This 
choice was made to properly normalize the Ha line pro- 
file for effective temperature determination (Sect. 13. 3p . 
Continuum normalization was done by fitting high or- 
der polynomials to the upper envelopes of each spectral 
order. Overlapping regions from continuous orders were 
later merged using the observed counts as weights. The 
order containing the Ha line as well as that containing 
the Mg I b lines were normalized differently. For these or- 
ders, we took advantage of the smooth variation of the 
blaze function from one spectral order to the next, to 
perform a 2D continuum normalization (i.e., both par- 
allel and perpendicular to the direction of dispersion) 
exclu ding orde r s affec ted by strong lines, as described in 
Bark lem et al.l ()2002l ). 

Daysky spectra were acquired in the afternoon of April 
24, 2011 using the same spectrograph but letting the 
sky light enter the slit from a "solar port" which points 
towards the zenith. The data reduction of these spectra, 
however, was different from that of the stars. Daysky fills 
up the slit as the source is not point-like. In addition, 
scattered light in the Earth's atmosphere is known to 
affect the line d epths and also equivalent widths (e.g., 
iGrav et all 120001 ). The latter effect is stronger for large 
angular separations between the Sun and the area of the 
sky observed and minimized when the area of the sky 
observed is very close to the location of the Sun. 
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Figure 1. Equivalent widths measured in our daysky spectra as 
a function of separation angle for three representative lines (open 
circles). Filled circles with error bars correspond to bin-averaged 
values. 

We obtained daysky spectra for zenith-Sun separation 
angles from 20 to 80 degrees. The measured equiva- 
lent widths clearly correlate with the separation angle 
(Fig. [1} while the angular variations seem more impor- 
tant for stronger features (Fig. ^. Our d aysky spec- 
tra re veal the systematic effects described in IGrav et aLl 
(|2000f ) and are thus not the best reference for very high 
precision abundance work. We obtained a number of 
daysky spectra when the Sun was at its highest point 
in the sky, about 20 degrees from zenith. Although not 
ideal, the problems described above are minimized under 
these circumstances. Thus, our solar reference spectrum 
corresponds to the average of 10 daysky spectra taken 
with the Sun near maximum height. We stress that the 
results obtained using our solar spectrum as reference are 
expected to be less precise than a differential analysis of 
16CygA relative to IGCygB, since the data acquisition 
and reduction process for the latter are more consistent. 

We can use the measurements from Fig. [T]to estimate 
the error in our measured line equivalent widths (EWs), 
since the S/N of each of our solar spectra is comparable 
to that of our 16Cyg A and B spectra. On average, the 
error is ~ 0.5 mA, which corresponds to about 1 % for a 
typical line used in this work {EW ~ 50 mA). 

We compiled a long list of spectral lines for the el- 
emental abundance determination. Due to the differ- 
ential nature of our work, and the similarity of the 
two brightest stars in the 16Cyg system, as well as 
their similarity to the Sun, accurate atomic data are 
not stri ctly necessary. We started with the linelist 
used bv lAsplund et al] (|2009l hereafter A09), which was 
carefully constructed for their solar abundance analysis. 
Not only are these lines mostly free from severe blends 
but also their atomic data are accurately known. We 
complemented the iron linelist with the compilation by 
[Ramirez et al.. (,20071 hereafter ROT), which consists of 
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Figure 2. Equivalent width difference measured in our solar spec- 
tra at 75 degrees minus 20 degrees of separation angle as a function 
of the equivalent width measured at 20 degrees of separation angle. 
Bin-averaged EWs (cf. Fig. [1} were used in this figure. 

lines that appear mostly unblended in solar-type stars 
and have reasonably accurate atomic data measured in 
the laboratory. For the rest of elements studied in this 
paper, we used additiona l lines listed in ot her large 
studie s of solar-type stars: iReddv et al.l (120(331. hereafter 
R03) iBensbv et all (f200l hereafter B05). l Neves"eraII 
(|2009l hereafter N09), and M09. Our adopted linelist 
is given in Table [H along with the source from which 
atomic data were adopted. When available, w e used 
the va n de r Waals damping constants by jBarklem et al.l 
(|200C1| ) and lBarklem fc Aspelund-Johanssoui (i2005i) . oth- 
erwise we adopted the values obtained from Unsold 's ap- 
proximation but with the Cg constants (e.g.. lGraylll992L 
page 217) multiplied by a factor of 3.1, a number that 
we estimated empirically from a study of the lines that 
are included in Barklem's work. 

Spectral line equivalent widths (EWs) were measured 
interactively using the splot task in IRAF. The locations 
of the local continuum, which are crucial for precise mea- 
surements of EW, were selected as carefully as possible. 
For each line, we searched for regions clearly free from 
spectral lines in the three spectra: solar, 16CygA and 
B. The exact same wavelengths for the continuum were 
used for all three spectra. If a small blend appeared to 
affect only one of the line wings, the other wing was used 
to measure EW by fitting only that side of the spectral 
line to a Gaussian profile. This procedure tends to re- 
sult in more precise EW values than deblending the line 
because of the uncertainty in determining the core wave- 
length of the blend. In these cases, we also ensured that 
the local continua for all three spectra were chosen as 
similarly as possible. We excluded features affected by 
telluric absorption. 

Equivalent width analysis allowed us to measure abun- 
dances of 22 elements: C, O, Na, Mg, Al, Si, S, K, Ca, 
Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, and Ba 
but we also measured abundances of the heavy elements 
Nd and Eu. Most of the Nd 11 and Eu 11 features are very 
weak and heavily blended in solar-type spectra. Thus, 
spectrum synthesis was necessary to measure their abun- 
dances. To ensure high precision, we restricted our anal- 
ysis to the best lines available, i.e., those with a very good 
continuum normalization and clearly identified blends 
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Table 1 

Adopted atomic data and our measured equivalent widths. 



A (A) 


species 


EP (eV) 


log{3/) 


source 


EW (mA) 
leCygA 


EW (mA) 
16 Cyg B 


EW (mA) 
Sun 


5052.1670 


C 


7.685 


-1.304 


A09 


40.0 


35.7 


34.7 


5380.3369 


C 


7.685 


-1.615 


A09 


27.2 


23.8 


22.7 


7771.9438 


O 


9.146 


0.352 


A09 


82.6 


74.8 


71.9 


7774.1611 


O 


9.146 


0.223 


A09 


71.5 


65.5 


62.5 


7775.3901 


o 


9.146 


0.002 


A09 


57.1 


51.7 


48.6 


6160.7471 


Na 


2.104 


-1.246 


A09 


61.2 


61.5 


55.9 


4730.0400 


Mg 


4.340 


-2.390 


R03 


75.9 


76.3 


69.1 


5711.0879 


Mg 


4.345 


-1.729 


A09 


111.6 


112.3 


107.8 


6318.7168 


Mg 


5.108 


-1.945 


M09 


45.2 


44.4 


43.7 


6319.2368 


Mg 


5.108 


-2.324 


A09 


31.6 


30.1 


25.0 


6696.0181 


Al 


3.143 


-1.481 


A09 


48.4 


47.9 


40.5 


6698.6670 


Al 


3.143 


-1.782 


A09 


25.9 


26.7 


21.2 


5665.5542 


Si 


4.920 


-1.940 


M09 


48.6 


47.7 


41.9 


5517.5400 


Si 


5.080 


-2.496 


N09 


17.6 


16.6 


13.7 


5690.4248 


Si 


4.930 


-1.773 


A09 


60.0 


58.0 


52.0 


5701.1050 


Si 


4.930 


-1.953 


A09 


46.0 


44.8 


38.8 


5772.1450 


Si 


5.082 


-1.653 


A09 


60.2 


58.1 


53.8 


5793.0732 


Si 


4.930 


-1.963 


A09 


52.2 


50.1 


45.6 


5948.5400 


Si 


5.080 


-1.208 


N09 


94.5 


92.6 


90.0 


6125.0298 


Si 


5.610 


-1.510 


R03 


39.3 


37.5 


32.8 


6142.4902 


Si 


5.620 


-1.540 


R03 


38.3 


37.5 


32.3 


6145.0200 


Si 


5.610 


-1.480 


R03 


47.8 


45.1 


39.6 


6237.3301 


Si 


5.610 


-1.116 


N09 


73.6 


71.4 


65.1 


6244.4800 


Si 


5.610 


-1.360 


R03 


55.7 


54.7 


48.1 


6721.8398 


Si 


5.863 


-1.060 


R03 


55.7 


54.0 


47.4 


7405.7700 


Si 


5.614 


-0.720 


B03 


101.4 


99.2 


94.8 


7800.0000 


Si 


6.180 


-0.710 


R03 


64.7 


60.1 


58.6 


4695.4429 


s 


6.525 


-1.829 


A09 


9.5 


8.2 


7.1 


6743.5400 


s 


7.866 


-0.600 


M09 


10.1 


8.4 


8.3 



with proper modeling and atomic parameters. We ana- 
lyzed 2 Euii lines (4205 and 6645. lOA) and 3 Ndii lines 
(5293.16, 5319.81, and 5811.57 A). For Eu we adopted the 
following isotopic fractions: ^^^Eu= 0.478 and ^^■^Eu^ 
0.522 (e.g., Lodders 2003). The relevant laboratory data 
used (transition probabilities, isot opic shifts, et c .) ar e 
given by Den Hartog et al. (200l. ILawler et all (|2001[ ). 
Ilvans et all (|2OO60 . and IRoederer et all (|2008D . Finally, 
lithium abundances were derived using line synthesis of 
the 67 08 A Li doublet, as described in iBaumann et all 

(poTol) . 

3. ATMOSPHERIC PARAMETERS 

The results presented in this section were ob- 
tained using the grid of MARCS model atmospheres 
(|Gustafsson et al.l 120081 ) available online^ Tests made 
with different flavors of Kurucz' model atmospheres show 
that the choice of atmospheric models has a very mi- 
nor impact on our derived relative parameters and abun- 
dances. We comment briefly on these effects in later sec- 
tions of this paper. Curve-of-growth analysis was used 
to measure elemental abundances from our equivalent 
width measurements using the 2010 version of the spec- 
trum synthesis code MOOG (Snedcn 1973)0 Line profile 
fitting was used to better constrain Tcs from the Ha line 
and logg from the Mgi6 lines. For the former, we used 
a fine grid of Ha theoretical line profiles (steps of 10 K 
in Toff, and 0.05 dex in logg and [Fe/H]), computed as 

^ http://marcs.astro.uu.se 

® http://www.as.utexas.edu/~chris/moog.htnil 



in iBarklem et al.l (|2002D Q Finally, the synthesis of the 
Mgi6 lines was done using SYNSPEC (Hubenv 1981 
iHubenv fc Lani 1995)0 

Stellar parameters often depend on each other so the 
process of their determination is intrinsically iterative. 
Here we present results which are self-consistent, in the 
sense that when we "adopt" a value below in the calcu- 
lation of a given parameter, we have made sure that no 
further changes to any of the other parameters involved 
will occur after another iteration. 

3.1. Excitation/ionization equilibrium 

The iron lines alone can be used to determine the 
atmospheric parameters using the conditions of ex- 
citation and ionization equilibrium. We start by 
assuming that the atmospheric parameters Teff (ef- 
fective temperature), logg (surface gravity), [Fe/H] 
(iron abundance) and Vt (microturbulent velocity) of 
16 Cyg A and B are solar (i.e., Tctf,logg,[Fc/B],vt = 
5777K, 4.437, 0.00, 1.00knis-i)I3 Then, Toff and Vt 
were iteratively modified so that trends of iron abun- 
dance (from Fei lines only) with excitation potential 



^ These Ha profiles were also computed using 
spheric models. 

* http://nova.astro.umd.edu/Synspec43/synspec.html 

^ We use the standard notation: Ax = log(nx/»iH) + 12, where 
rix is the number density of element X, and [X/H] = Ax — A^. 

^^ The solar microturbulent velocity vt is somewhat uncertain 
and dependent on the iron line list used and model atmosphere 
adopted. Its exact value is not crucial for our strictly differential 
study, as shown later. 
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Figure 3. Top panels: [Fe/H] values deterinined on a linc-by-line 
basis witli respect to our solar spectrum as a function of excitation 
potential for 16 Cyg A and B; Fel (Fen) lines are shown with open 
(filled) circles. The dot-dashed line is at [Fe/H] = +0.104 while 
the dashed line is at [Fe/H] = +0.057, which are the mean [Fe/H] 
values for 16 Cyg A and B, respectively. Bottom panels: as in the 
top panels but as a function of reduced equivalent width. 

(EP) and reduced equivalent width {log EW/X), respec- 
tively, approached zero. Simultaneously, log 5 was mod- 
ified so that the mean difference in iron abundance from 
Fel and Fell lines also approached zero. Errors in the 
derived parameters were computed from the uncertainty 
in the measured slope of abundance versus EP (for T^e) 
and from the standard deviation of the mean difference 
of Fei and Fen abundances (for log 5). 

Using the excitation/ionization equilibrium method we 
obtain: T^s = 5818 ± 38 K, logg = 4.31 ±0.05, [Fe/H] = 
-h0.104 ± 0.018 for 16 Cyg A; and Tch = 5742 ± 30 K, 
log(7 = 4.31 ±0.05, [Fe/H] = +0.057 ±0.017 for 16 Cyg B 
(see Fig. [3]) . If we use Kurucz no-overshoot models in- 
stead of MARCS we find a shift of only +5K in Toff, 
+0.01 dex in \ogg, and +0.002dex in [Fe/H], system- 
atically consistent for both stars. JVIoreover, adopting 
a slightly different solar microturbulence: vt = 1.0 ± 
O.lkms"^, we find changes of only ±0.003dex in [Fe/H] 
and, as expected, shifts of ±0.1 km s""'^ in the microtur- 
bulent velocities of 16 Cyg A and B. This shows that the 
choices of model atmosphere flavor and solar microturbu- 
lence have an almost negligible impact on our differential 
abundance results. 

The calculation described in this section suggests that 



16 Cyg A is slightly more metal-rich than 16 Cyg B. How- 
ever, since there is a non- negligible difference in the ef- 
fective temperature of the stars, small differential sys- 
tematic effects might still be present and it is therefore 
important to check the parameters here derived using 
other independent methods. 

3.2. Photometric effective temperature 

ICasagrande et al.l (J2010D have implemented the so- 
called infrared flux method (IRFM) to derive the effec- 
tive temperatures of a large sample of dwarf and sub- 
giant stars. The zero point of their temperature scale 
has been determined with an unprecedented precision 
of only 15 K. Metallicity-dependent color-Toff transfor- 
mations provided in that work allow us to estimate Toff 
values for the 16 Cyg A and B pair from their observed 
photometry. We used the {B — V) and {b — y) colors 
and our derived [Fe/H] values to obtain: Toff = 5792 K 
for 16 Cyg A and Toff = 5733 K for 16 Cyg B. Photomet- 
ric errors translate into an error of about 16 K for both 
stars, similar to the color-to-color error. The uncertain- 
ties in [Fe/H] translate into a 5K uncertainty in Teff. 
Adding these errors in quadrature we estimate a random 
error of 23 K, which has to be added l inearly to the sys- 
tematic error in the zero point of the ICasagrande et ahl 
(pOlO) Toff scale to estimate the uncertainties in our de- 
rived Toff values. Thus, the effective temperatures we 
derived using colors have errors of 38K. 

The color temperatures, which are based on the most 
reliable IRFM temperature scale to date, confirm that 
16 Cyg A is warmer than 16 Cyg B, although the differ- 
ence seems to be smaller when using the IRFM Toff scale 
(59 K) instead of the spectroscopic method described in 
the previous section (76 K). Nevertheless, the Teff values 
from the two methods are in good agreement within the 
la errors. 

3.3. Model fits to the Ha line profile 

The wings of the Balmer lines are known to be sen- 
sitive to Teff and if \ogg and [Fe/H] can be estimated 
with high precision independently, as in our case, model 
fits to the Ha line profile can provide very precise Tpft 
values (e.g..lFuhrmann et al.|[T99l IBarklem et al.ll20C)2l: 
iRamirez et al.ll2006D . This is particularly the case if a 
good continuum normalization is made using very high 
quality spectra. 

We used least squares minimization to find the best 
model fits to our Ha data as well as the formal errors 
in Toff. Careful inspection of our spectra allowed us to 
find several wavelength windows free from spectral lines, 
both stellar and telluric, other than Ha, which were then 
used in the model fits. The same windows were used for 
the three spectra. In Fig. |4] we show the observed nor- 
malized Ha lines for 16 Cyg A, B, and the Sun (daysky) 
along with their best fit models. The problems described 
in Sect. [5] for the solar spectrum are not as important for 
this broad feature, in particular because we are only in- 
terested in the wings, so we can use reliably our daysky 
spectrum as solar reference in this case. For our solar 
data, we derive Teff = 5732 ± 32 K, similar to the value 
obtained bv IBarklem et al.l ()2002| )R This suggests that 

^^ Our error bar here corresponds to observational noise only. 
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Figure 4. Observed normalized Ha line profiles in our 16CygA, 
16CygB, and daysky (Sun) spectra. Best fit models for these lines 
are also shown. 

there is still roora for improvement in the modeling of 
Balmer lines. Hereafter, we apply a zero point correction 
of +45 K to the effective temperatures derived from Ha 
since this brings our derived solar Tcff up to the nominal 
value of 5777 K. 

We computed T^s values from the Ha line profile fits 
for a small grid of \ogg = 4.2,4.3,4.4 and [Fe/H] = 
0.00,0.05,0.10. Then, we used the grid of Tcff results to 
interpolate to our preferred log g and [Fe/H] parameters 
as well as to compute the error introduced by uncertain- 
ties in log (7 and [Fe/H]. With our final \ogg and [Fe/H] 
values we obtain: T^s — 5819 ± 25K for 16CygA and 
Teff = 5762 ± 26K for 16CygB. An error of 0.02 dex in 
log g translates into only 2 K whereas an error of 0.03 dex 
in [Fe/H] translates into about 1 K of Toff. The uncer- 
tainties in our derived values of logg and [Fe/H] thus 
have negligible impact on this calculation. Moreover, 
given the very high signal-to-noise ratios of our spectra 
in this region {S/N ~ 500, obtained because the spec- 
trograph setup was chosen to maximize S/N at 6562 A), 
our Toff values derived from Ha line profile fitting are the 
most reliable, as the smaller errors in comparison with 
Toff values derived with other methods suggest. 

Our three estimates of Toff arc in excellent agreement 
with each other for both 16CygA and B. Since these 
techniques are independent, we can obtain a final Toff 
value by averaging them, using the estimated errors as 
weights. We find: Teff = 5813 ± 18K for 16CygA and 
Teff = 5749±17K for 16CygB0The degree of precision 
with which we have been able to measure the effective 
temperatures of these stars is almost unprecedented for 



IBarklem et al.l II2002I ) studied systematic errors very carefully and 
obtained an error of about 80 K for their solar analysis. To remain 
consistent with the rest of our strictly differential work, we consider 
the observational errors only. 

^^ We calculate the weighted mean value {x) and error (a) of N 
data points Xi with errors ai as follows: x = '^{xi/cr^)/ ^{1/cr^), 
a^ = 1/ 5I](1/'''?)- We note, however, that this formula for a does 
not take into account possible systematic errors and could also 
under- or over-estimate the error if the ctj values are not realistic. 
The error computed from the sample variance is in principle more 
reliable: u"^ = {E[(a:» - 2;)V<^?)]/ E(l/^?)}/{A' - !)■ Using the 
latter we find Teff errors of 8 and 9K for i6CygA and B, respec- 
tively, but we adopt the former, more conservative error bars. 



any other star but the Sun. The only other exception 
is the case of the sol ar twin star HIP 56948, for which 
iMelendez et al.l (|2011f) report a Toff with an error of 7K. 

3.4. Trigonometric surface gravity 

The parallax of the 16 Cyg pair has been accurately 
meas ured by the Hipparcos mission (iPerrvman et aLl 
|1997( ). The new reduction bv Ivan LeeuwenI (|2007D gives 
values of 47.44 ± 0.27mas for 16CygA and 47.14 ± 
0.27 mas for 16 Cyg B. The V magnitudes of these stars 
are also precisely known: V{A) = 5.959 ±0.009, V(B) = 
6.228±0.019. Using these numbers we can infer the abso- 
lute V magnitudes of the stars: My (A) = 4.340 ±0.015, 
My(B) = 4.595 ± 0.023. Since their effective tempera- 
tures are also known with very high precision, we can 
reliably place the stars on the HR diagram and compare 
their location with theoretical stellar evolution predic- 
tions (isochrones) to infer their masses, ages, and surface 
gravities. 

Our particular implementation of stellar parame- 
ter determinat i on us ing isochron es is described i n 
IMelendez et al.l ()201lL see also Ba umann et al.l |2010[ ) . 
Our method follows the probabilistic scheme explained in 
R03 (their Sect. 3.4) and also All ende Pri cto et al. C200j, 
their Appendix A) but the choice of isochrones is differ- 
ent. Basically, all isochrone points within a radius of 
3(7 from the observed parameters are used to compute 
mass, age, and logg probability distribution functions, 
from which a most likely value and 68% (Icr-like) and 
95% (2cr-like) confidence limits ca n be derived. Finely 
spaced Yonsei-Y ale isochrones (e.g.. lYi et al.ll200Tl 120031 : 
iKim et 311120021) were used in these computations. 

Using our isochrone method we infer masses of Af (A) = 
1.05 ± 0.02 Mq and M(B) = 1.00 ± 0.01 A/©. For the 
surface gravities we find: \ogg{A) = 4.28 ± 0.02 and 
log (7(B) = 4.33 ± 0.02. Since the two stars were born 
at the same time, it is expected for the primary to be 
slightly more evolved given its somewhat larger mass. 

The ages that we infer for the stars in the 16 Cyg sys- 
tem are: t(A) = 7.151^3 Gyr and t(B) = 7.26+;^:^!^ Gyr. 
Error bars correspond to the Icr-like lower and upper lim- 
its. The age probability distributions (APDs) for these 
stars are shown in Fig.[S] The asymmetry of these APDs 
explains the very different upper and lower la errors 
given above. Although they are somewhat broad, the 
APDs of 16 Cyg A and B peak at nearly the same age. 
A much more precise age can be obtained by combining 
the two APDs, as shown with the solid line histogram in 
Fig. m This combined 16 Cyg APD is simply the prod- 
uct of the APDs of the components of the 16 Cyg binary 
system, normalized to 1 at the maximum. Using the com- 
bined APD we infer an age of 7.10to'35Gyr for 16 Cyg. 
Note that this very small age error bar does not include 
systematic errors due to model uncertainties or errors in- 
troduced by statistical biases (see below). It represents 
only the internal precision of our method when applied 
to a binary system. 

Stellar parameters derived from isochrones are sub- 
ject to a number of systematic errors as well as sta- 
tistical biases due to isochrone sam pling which are 
descri bed in detail by, for example, iNordstrom et al.l 
(|2004 their Sect. 4.5.4). Bayesian approache s have been 
proposed to tackle these problems (e.g., ,Pont fc Eveii 
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Figure 5. Normalized age probability distributions for 16 Cyg A 
(dashed histogram) , 16 Cyg B (dot-dashed histogram) , and the 
16 Cyg binary system (solid line histogram) . 

20P; ■j0rEensen & LindegrenllMl Ida Silva et all [20061: 
Casagrande et al. 2011). We can comp are our results 



with the Bayesian implementation by Ida Silva et al.l 
(|2006f ) since they provide a very convenient online tool 
that uses their methodEB We find log 5(A) = 4.26±0.0f 
and log (7(B) = 4.3f ± O.Of , in excellent agreement with 
our estim ates even thou gh the choice of isochrones is 
different (Ida Silva et al.l use PADOVA isochrones; e.g., 
IBertelh et al.l 119941 : IGirardi et al.ll2000D and we do not 
cor rect for statist i cal bi ases. The masses derived with 
the Ida Silva et al.l ()2006t ) online tool are about 0.02 Mq 
lower but the difference betwee n 16 Cyg A and B is nearly 
identical. For the stellar ages, Ida Silva et al.l online tool 
suggests 6.7 ± 0.5 Gyr (A) and 7.6 ± 0.7 Gyr (B). These 
values are also consistent with our derived age for the 
16 Cyg system considering the la uncertainties. 

Note that the excitation/ionization method of atmo- 
spheric parameter determination was not able to resolve 
the small difference in log g found in this section although 
the results are also consistent within Icr. Clearly, for 
nearby stars with well determined effective temperatures, 
parallaxes, and good photometry, the isochrone method 
is superior. 

3.5. Surface gravity from the Mgih lines 

The wings of certain strong spectral features are mildly 
sensitive to the stellar surface gravity. A particularly 
important feature in this regard is the Mg i h triplet at 
about 5180 A. The effective temperature and magnesium 
abundance must be known with high precision in order 
for this line to be useful in the determination of log (7. 

Fig.|6]shows the spectral region around the strongest of 
the Mg I h lines and the red wing of the middle one along 
with synthetic spectra computed for values of log g that 
reproduce best the regions free from lines other than the 
Mgi6 lines, such as those around 5174 A and 5185 A. We 
fine-tuned the transition probabilities of these lines to 
reproduce the solar Mg i h lin es with a Mg abunda nce of 
AMg = 7.60, as measured bv lAsplund et al.l ()2009[) . For 
16 Cyg A and B we used synthetic models with our final 
Toff and [Fe/H] values as well as our derived magnesium 
abundances. Synthetic profiles for logg±0.2dex are also 

^'^ http;//stev. oapd.inaf.it/cgi-bin/param 
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Figure 6. Spectral region showing the strongest of the Mg I b lines 
and the red wing of the middle Mglfe line in 16 Cyg A, 16 Cyg B, 
and the Sun. Dashed lines show best fit models for each star (with 
their derived logg values) while dotted lines correspond to syn- 
thetic spectra with logg ±0.2. 



Table 2 

Effective temperature and surface gravity determinations 



Teff (K) 


16 Cyg A 


16 Cyg B 


Fe lines 
Colors (IRFM) 


5818 ± 38 
5792 ± 38 

5819 ± 25 


5742 ± 30 
5733 ± 38 
5762 ± 26 


Adopted 


5813 ± 18 


5749 ± 17 


logg [cgs] 


16 Cyg A 


16 Cyg B 


Fe lines 

Trigonometric 

Mgib 


4.31 ±0.05 
4.28 ±0.02 
4.27 ±0.04 


4.31 ±0.05 
4.33 ±0.02 
4.33 ±0.04 


Adopted 


4.282 ±0.017 


4.328 ±0.017 



Table 3 

Basic data and derived parameters for the 16 Cyg binary system 





16 Cyg A 


16 Cyg B 


V (mag) 


5.959 ± 0.009 


6.228 ±0.019 


vr (mas) 


47.44 ± 0.27 


47.14 ±0.27 


Mv (mag) 


4.340 ±0.015 


4.595 ± 0.023 


Tefr (K) 


5813 ± 18 


5749 ± 17 


logg [cgs] 


4.282 ±0.017 


4.328 ±0.017 


Mass (Mq) 


1.05 ±0.02 


1.00 ±0.01 



shown in Fig. [6] to illustrate the sensitivity of these fea- 
tures to logg. From a least-squares minimization of the 
clean regions we find logg — 4.27 ±0.04 for 16 Cyg A and 
logg = 4.33 ± 0.04 for 16 Cyg B, in excellent agreement 
with the values found in the previous section. 

Similar to the case of Toff, the three methods used 
to derive logg are independent so a weighted average 
can be used to determine the best solution. We find: 
logg(A) = 4.282±0.017, logg(B) = 4.328±0.017. These 
and other fundamental parameters derived in this section 
are summarized in Tables [2] and |3l 

4. ELEMENTAL ABUNDANCES 

Curve of growth analysis within MOOG was used to 
determine elemental abundances from equivalent width 
measurements of lines in the spectra of 16 Cyg A, B, and 
the Sun. We measured 288 lines for 22 elements: their 
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1500 



Figure 7. Elemental abundances of 16CygA and B relative to 
our measured solar abundances as a function of dust condensation 
temperature for light elements (Z < 30). Solid lines show the 
mean trend of solar twin stars according to Melcndez ct al. ( 2009) 
but shifted upwards in each case to match the mean abundance of 
volatile {Tq < 900 K) elements. Open circles correspond to Mn, 
Co, and Sc II. 

EWs are given in Table [H Since the oxygen abundances 
inferred from the infrared triplet lines, which are the 
only oxygen lines we u se, arc severely affected by non- 
LTE effects (e.g., rKiselman 1993; Fabbian ct al. '2009), 
and reliable non-LTE computations can be found in the 
literature, we used the non-LTE correction tables by 
iRamirez et al.l ()2007[ ) to take this systematic effect into 
account. Lines from two species (neutral and singly ion- 
ized) of Sc, Ti, Cr, and Fe were available; for the other el- 
ements only one species was measured. Good agreement 
was found for the abundances derived from two species of 
the same element except for Sc (see below). The abun- 
dances of Nd and Eu were determined using spectrum 
synthesis of a few carefully selected features. Spectrum 
synthesis was also employed to measure lithium abun- 
dances from the 6708 A Li doublet. 

We measured differential abundances on a line-by-line 
basis, which reduces the impact of errors in the atomic 
data and systematics in the modeling, given that these 
are all similar stars. Abundance ratios relative to solar 
([X/H]) are listed in Table S) Errors in this table cor- 
respond to the standard error of the line-to-line scatter 
added in quadrature with the random error introduced 
by the uncertainties in our derived atmospheric parame- 
ters. Fig. [7] shows the [X/H] abundances as a function of 
the du st condensation temperature sequence by Lodders 
(|2003[ ) for the light elements (Z < 30). Heavy element 
{Z > 30) abundances are shown separately in Fig. [5] In 
both Figs. [7] and [5] we plot with a solid line the mean 
trend of solar twin stars by M09, shifted to match ap- 
proximately the abundance pattern of Z < 30 elements 
for each star. 

The weighted mean [X/H] abundances from all species 
studied, which we denote by [M/H], are -h0.103 ± 0.023 
for leCygA and -F0.069±0.026for 16CygB. The 16Cyg 
binary system is therefore more metal-rich than the Sun 
but, interestingly, a small difference in the overall metal 
content appears to make 16CygA more metal-rich than 
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Figure 8. As in Fig.[7]for the heavy elements [Z > 30). Stars (cir- 
cles) denote elemental abundances measured with spectrum syn- 
thesis (equivalent widths). 



16 Cyg B. 

Considering only the light elements {Z < 30, Fig. [7]), 
the observed abundances of 16 Cyg A and B follow the 
typical pattern of solar twin stars, i.e., they correlate 
with dust condensation temperature as in M09. Note 
that the latter was obtained from the average abundances 
of 11 solar twin stars, resulting in a very small abundance 
error (of order 0.01 dex). However, for individual stars, 
the dispersion seen in Fig. [7] is fully consistent with the 
analysis by M09. 

On the other hand, the heavy element abundances 
{Z > 30, Fig. HI) do not follow the typical solar twin 
behavior. Strictly speaking, 16 Cyg A and B are solar 
analogs and not solar twins, mainly due to their slightly 
lower surface gravities, which makes them more evolved 
than the Sun. Compared to the Sun, the 16 Cyg sys- 
tem is about 2.5 Gyr older and it is therefore likely that 
the anomalous abundances of the heavy elements (when 
compared to solar twins) are due to Galactic chemical 
evolution effects. We discuss these in detail in Sect. 15.21 

In principle, the mean metallicities ([M/H]) of 
16 Cyg A and B are consistent within Icr. We should 
point out, however, that this analysis uses our solar spec- 
trum as reference. Although obtained with the same 
instrument, which is ideal, our solar spectrum was not 
observed in the exact same way the stellar spectra were 
taken. Important differences in the data reduction proce- 
dures are also present. Furthermore, our reference spec- 
trum is based on observations of the daysky, which is 
known not to be an exact replica of the solar spectrum, 
although we attempted to minimize its impact by using 
daysky spectra at small angular distances from the Sun, 
as explained in Sect. [5] 

Inspection of Figs. [7] and [8] shows clearly that the ele- 
ment to element scatter is very similar in 16 Cyg A and B. 
This suggests that the difference between A and B might 
be better established by comparing them directly rather 
than using our possibly faulty solar spectrum. Galac- 
tic chemical evolution effects will also be minimized, if 
not eliminated, with this approach. Thus, we measured 
differential abundances of 16 Cyg A, again on a line-by- 
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Table 4 

Elemental abundances relative to solar ([X/H]) and differential between 16CygA and B {A[X/H]). 





Species 


Tc (K) 


[X/H] 


error 


[X/H] 


error 


A [X/H] 


error 








(leCygA) 


(leCygB) 


(leCygA- 


-leCygB) 


6.0 


C 


40 


0.043 


0.012 


0.007 


0.012 


0.037 


0.016 


8.0 


O 


180 


0.075 


0.016 


0.052 


0.018 


0.024 


0.024 


11.0 


Na 


958 


0.102 


0.032 


0.068 


0.031 


0.034 


0.020 


12.0 


Mg 


1336 


0.104 


0.026 


0.060 


0.026 


0.043 


0.013 


13.0 


Al 


1653 


0.138 


0.018 


0.110 


0.018 


0.028 


0.023 


14.0 


Si 


1310 


0.118 


0.006 


0.074 


0.004 


0.044 


0.008 


16.0 


S 


664 


0.064 


0.029 


0.020 


0.030 


0.044 


0.016 


19.0 


K 


1006 


0.124 


0.034 


0.076 


0.034 


0.048 


0.028 


20.0 


Ca 


1517 


0.118 


0.019 


0.074 


0.018 


0.044 


0.017 


21.0 


Scl 


1659 


0.115 


0.065 


0.041 


0.065 


0.074 


0.027 


21.1 


Sell 


1659 


0.139 


0.013 


0.098 


0.012 


0.041 


0.011 


22.0 


Til 


1582 


0.109 


0.018 


0.076 


0.016 


0.033 


0.023 


22.1 


Till 


1582 


0.124 


0.017 


0.084 


0.017 


0.040 


0.009 


23.0 


V 


1429 


0.120 


0.020 


0.083 


0.020 


0.037 


0.026 


24.0 


CrI 


1296 


0.085 


0.014 


0.053 


0.013 


0.033 


0.018 


24.1 


CrII 


1296 


0.095 


0.010 


0.061 


0.011 


0.034 


0.013 


25.0 


Mn 


1158 


0.107 


0.029 


0.093 


0.028 


0.014 


0.022 


26.0 


Fel 


1334 


0.104 


0.012 


0.061 


0.011 


0.042 


0.016 


26.1 


Fell 


1334 


0.089 


0.009 


0.064 


0.010 


0.025 


0.013 


27.0 


Co 


1352 


0.100 


0.018 


0.093 


0.016 


0.007 


0.018 


28.0 


Ni 


1353 


0.109 


0.012 


0.073 


0.010 


0.036 


0.014 


29.0 


Cu 


1037 


0.082 


0.028 


0.032 


0.028 


0.051 


0.019 


30.0 


Zn 


726 


0.120 


0.033 


0.066 


0.033 


0.053 


0.008 


39.1 


Y 


1659 


-0.012 


0.033 


-0.048 


0.044 


0.036 


0.016 


40.1 


Zr 


1741 


0.061 


0.031 


0.022 


0.031 


0.039 


0.018 


56.1 


Ba 


1455 


0.034 


0.022 


-0.006 


0.039 


0.040 


0.020 


60.1 


Nd 


1594 


0.080 


0.021 


0.043 


0.032 


0.030 


0.035 


63.1 


Eu 


1347 


0.125 


0.035 


0.095 


0.035 


0.037 


0.023 



line basis, but using 16CygB as the reference star. We 
denote these differential abundances with the symbol 
A [X/H]. The results are given in Table ID Similar to the 
errors in [X/H], the errors listed in Table |4] for A [X/H] 
correspond to the standard error of the line-to-line scat- 
ter added in quadrature with the errors introduced by un- 
certainties in our derived stellar parameters. Fig.[S]shows 
the derived 16CygA - 16CygB relative abundances as 
a function of atomic number and dust condensation tem- 
perature. Note that in this case both the light and heavy 
elements behave similarly. 

Fig. [9] shows that not only are the la errors of the 
A-B relative abundances small but it also reveals that 
the element-to-element dispersion is extremely small. 
A weighted average of all species measured suggests 
a metalhcity offset of A[M/H] = -t-0.040 ± 0.010 dex 
for 16CygA relative to 16CygB. The three species 
that depart the most from this average value are Sc i, 
Mn, and Co. It is tempting to attribute these dis- 
crep ant data poin ts to strong differential non-LTE effects 
(see iZhang et al.l [2008 for Sc i . Bergemann fc Gehreiil 
20071 for Mn. and~lBergemannl 120081 for Co: see also 
Bergemann et all 120101 ) since there is a difference of 
about 65 K in Tes between the two stars in addition to a 
difference of 0.05 dex in log 5. Furthermore, we note that 
the Sen abundance, which is known to be more reliable 
than Sci, is consistent with the derived A[M/H] while 
the Icr line-to-line scatter (not the standard error) of the 
Mn abundance difference is the largest of all species an- 
alyzed with equivalent widths. A more proper treatment 
of Co line formation requires consideration of hyperfine 
structure (HFS) splitting, which we tested for three Coi 



features using HFS constants by iKurucz fc: Belli ()1995D . 
In this case we obtained A[Co/H] = -hO. 007 ± 0.018, in 
excellent agreement with the value derived from equiv- 
alent widths and ignoring HFS. Perhaps a combination 
of HFS and non-LTE effects cou ld explain our Co abun- 
dances. iBergemann et all (|2010l ). for example, calculate 
a non-LTE correction of 0.13 dex for the Sun, which is 
relatively large. Excluding arbitrarily these three dis- 
crepant data points (Sci, Mn, and Co), which are likely 
affected by small but non-negligible systematic errors in 
high precision work, the metalhcity offset is: A[M/H] = 
-1-0.041 ± 0.007 (A-B). Hereafter, we discuss our results 
ignoring the Sc i, Mn, and Co abundances. 

Note that the la uncertainty of the metalhcity differ- 
ence that we inferred above corresponds to about 1 % on 
a linear scale. Interestingly, the Icr error of our mea- 
sured EWs is also about 1 %. This suggests that the 
error in our very precise elemental abundance determi- 
nations, which are based on curve of growth analysis (ex- 
clusively on the linear part) with the EWs as input data, 
is dominated by the observational errors (elements Nd 
and Eu are the only exceptions, since their abundances 
were measured using line profile fitting). Thus, it should 
be in principle possible to improve upon these estimates 
using even higher quality data. 

We have explored the impact of systematic errors in 
our derived atmospheric parameters by arbitrarily mod- 
ifying them and recomputing the abundances. First, 
we tried increasing and decreasing the Tog values of 
both stars by 100 K. Interestingly, in both cases we find 
that the A-B metalhcity difference decreases to about 
0.03 dex but the element-to-element scatter increases to 
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Figure 9. Top panel: elemental abundance difference between 
16 Cyg A and B as a function of atomic number. Open symbols 
show the three species more discrepant from the mean: Scl (21), 
Mn (25), and Co (27). Bottom panel: as in the top panel for 
the abundance differences versus dust condensation temperature. 
The dashed line is at +0.041 while the solid lines represent the 
mean trend of solar twins by M elcndez et al.l 1)20091 ) with two ar- 
bitrary offsets. The dot-dashed line corresponds t o a s lope of 
1.4 X 10-5 dexK^i, as derived bv lLaws fc Gonzaled |[200!D . 



0.008 dex. Adopting identical log g = 4.30 values for both 
stars increases the metallicity difference to -1-0.046 dex 
but in this case the element-to-element scatter increases 
significantly to 0.013 dex. The fact that the element- 
to-element scatter naturally minimizes for our derived 
parameters suggests that they are very reliable. 

We have also re-determined our ab undance s usin g 
the differential parameters derived by iTakedal ()2005[ ). 
who claims that there is no metallicity difference be- 
tween 16 Cyg A and B. We kept our 16 Cyg B parameters 
fixed an d adopted 16 Cyg A parameters consistent with 
iTakedaf s. The inferred metallicity difference in this case 
is -h0.032 ± 0.013, i.e., we find that 16 Cyg A is still more 
metal-rich than 16 Cyg B. However, for these parame- 
ters we detect a small negative trend for the abundance 
differences and condensation temperature which makes 
the difference smaller for refractory elements like iron. 
For those elements the metallicity difference is about 
-1-0.02 ±0.01 dex. It seems possible that a combin ation of 
these effects is responsible for the ITakedal (j2005[ ) results 
but we admit that we are unable to reproduce exactly 
his derived elemental abundances. 

Our main result that 16 Cyg A is more metal-rich than 
16 Cyg B was already suggested in the high precision 
work by iLaws fc GonzalezI ()200lD . They found, how- 
ever, a metallicity difference of -1-0.025 ± 0.009 based 
only on the analysis of iron lines. As a note added in 
proof they commented that abundances of 13 other ele- 



ments were measured and showed a small positive cor- 
relation with dust condensation temperature. Their de- 
rived slope is shown in the bottom panel of our Fig. |9l 
Depending on the elements that one uses to calculate 
the slope, it could be poss ible to find goo d agreement 
between our data and the iLaws fc Gonzaled result. Nev- 
ertheless, using our abundances we obtain essentially 
a zero slope (—0.2 ± 0.6 x 10^ dex K "^), which is in- 
consistent with the slope derived by ILaws fc GonzalezI 
(+1.4 ± 0.5 X lO^dexK"-^). We argue that our results 
are more robust owing to the larger number of elements 
analyzed and the very high precision achieved for the 
overall metallicity difference. 

Finally, using Kurucz' model atmospheres we find 
A[M/H] = +0.039 ± 0.007 (no-overshoot models) and 
A[M/H] = +0.040 ±0.007 (overshoot models). Thus, we 
demonstrate again that the choice of model atmosphere 
grid has an almost negligible impact on our results. 

5. DISCUSSION 
5.1. Planet signatures in stellar abundances 

High precision elemental abundance work in solar twins 
and analogs has revealed a possible connection between 
the properties of a star and the process of planet forma- 
tion that occurred around it. The work by M09 showed 
that, compared to the majority (~85%) of solar twin 
stars, the Sun exhibits a deficiency of refractory elements 
relative to volatiles (of order 20%). They proposed that 
this solar abundance anomaly could be due to the for- 
mation of the terrestrial planets in the solar system since 
these objects are rich in refractory elements. Thus, the 
missing mass of refractory elements in the solar photo- 
sphere would be currently located in the terrestrial plan- 
ets. This hypothesis assumes that the majority of so- 
lar twin stars did not form as many rocky objects as 
in the solar system or that if they did, they were ac- 
creted by the star at later times, possibly through orbital 
migration. Only about 15 % of solar twin stars show a 
photospheric composition indistinguishable from the so- 
lar one. These objects are potentiall y terrestrial p lanet 
hosts. Later work by R09 and Ramirez et al.l (J201C() con- 
firmed the observational results of M09 and expanded the 
sample towards higher metallicities, where they found 
that, based on the abundances measured, the fraction of 
terrestrial planet hosts increases towards higher [Fe/H] . 

Thus, it seems likely that a "planet signature" has been 
imprinted in the chemical composition of the host stars 
and, arguably, it appears to relate only to the formation 
of rocky objects. Indeed, M09 showed that metal-rich so- 
lar analogs with giant planets detected have abundance 
patterns relative to iron similar to those of the major- 
ity of solar twins, suggesting that the solar abundance 
anomaly could not be related to the formation of the gas 
giants. This type of analysis, however, should be exam- 
ined more carefully because the impact of systematic er- 
rors and Galactic chemical evolutio n could be importan t 
in this higher [Fe/H] regime (see also lSchuler et al1l2011|) . 
Data by independent groups have confirmed the obser- 
vational result by M09 and R09 ([R amirez et al. 201^ 
iGonzalez Hernandez et al1l2010HGonzalez et al..,2010bi) . 
Although the planet forma tion signature interpreta- 
tion h as been contended by IGonzalez Hernandez et ahl 
(|201Clf ). a response to their criticism has been provided by 
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iRamirez at alj (|2010[ ). Alternative scenarios to explain 
the solar abunda nce anomaly have als o been explored, 
without success (IR amirez et al.l l2010t iMelendez et alJ 
120111 : iKiselmaneta l. 2011). 

The planet signature hypothesis described above has 
been quantitatively tested by Chambers (2010), who 
computed the effect of adding two Earth masses (2M0) of 
Earth-like material and 2Mq of CM chondrite material 
to the present-day solar convective envelope. He finds 
that in this experiment the solar abundance anomaly 
would disappear, i.e., the Sun would not show the de- 
ficiency of refractory elements anymore. Thus, his cal- 
culations support the hypothesis that the missing mass 
of rock in the solar photosphere is currently locked-up in 
the rocky bodies of the solar system. A simi l ar est imate 
has been recently made by Mclcnde z et al.l (|2011f ) who 
suggest an even larger amount of CM chondrite mate- 
rial (4M0). Although such large mass is not observed 
in the present-day asteroid belt, it has been suggested 
that small rocky bodies were very abundant when the 
Sun was formed, but most of the material from the as- 
teroid belt was probably removed by strong resonances 
(e.g.,|Wcidcnschi lhng 1977; Chambers & WcthcriU 2001; 
iPetit et al.ll200l . 

Excluding the heavy elements (i.e., for Z < 30), the 
abundance patterns of both 16CygA and 16CygB are 
compatible with those observed in most solar twin stars 
(Fig. |7|). Within the planet signature scenario, this im- 
plies that none of these objects host terrestrial planets 
or that any rocky body that once formed there has been 
accreted by the host star. Perhaps binary systems like 
16 Cyg are not stable enough to retain any terrestrial 
planets that might form. The gas giant planet around 
16 Cyg B is highly eccentric, a property that has been 
attributed by some a uthors to the influen c e of the stel- 
lar c ompanion (e.g., iCochran et al.l ITQQTI : iMazeh et al.l 
Il997f ). The lack of rocky planets around 16 Cyg B could 
also be related to the presence of the gas giant planet. 
As evidenced by the Kepler mission data, multiple sys- 
tems with small pl anets rarely contain hot Jupiters 
([Latham et al.ll2011| ). In any case, the abundances of 
16 Cyg A and B are not unusual considering that only 
about 15 % of solar twins share the exact same (solar) 
composition. 

5.2. Impact of Galactic chemical evolution 

The heavy elements clearly do not follow the typical 
solar twin trend (Fig. |8|). With the exception of Eu, 
the other Z > 30 elements are found significantly be- 
low the mean solar twin trend. The Nd abundances are 
marginally consistent with the expected trend but the Y, 
Zr, and Ba abundances are irreconcilable with that pat- 
tern. These elements are produced primarily by neutron- 
capture reactions on timescales slow (the s-process) or 
rapid (the r-process) relative to the average /9-decay rates 
of radioactive isotopes along the reaction chain. The s- 
process is associated with low- and intermediate-mass 
stars (1.3 ^ Af < 8.0 Mq) on the asymptotic giant 
branch. The site (or sites) of the r-process is not known 
but is suspected to be associated with core collapse su- 
pernovae (M > 8.0 Mq). These stellar timescales imply 
that r-process enrichment occurs earlier than s-process 
enrichment, which is confirmed by the observed heavy 
element abundance patterns in low metallicity stars in 



the Galactic halo and disk (e.g.. lMcWilliamlll997l) . 

Each of the r- and s-process contribute about half of 
the elements with .Z > 30 i n the S un, and, following pi- 
oneering work bv lCameronI (|1973l ). the solar abundance 
of each isotope (or element) can be decomposed into r- 
and s-process fractions. For example, 72% of the solar 
Y wa s produced via the s-process (e.g., iSimmerer et al.l 
[2OOI . Similarly, 81% of Zr, 85% of Ba, 58% of Nd, and 
3% of Eu in the Sun were produced by the s-process. 
These fractions represent the chemical inventory of one 
point in the Galactic interstellar medium 4.5 Gyr ago. 
The stars in the 16 Cyg system formed 2.5 Gyr earlier 
and they are slightly deficient in s-process material rela- 
tive to the Sun. As can be seen in Fig. |8l elements with 
a significant s-process component, such as Y, Zr, and 
Ba, are the most deficient. Nd, which has roughly equal 
contributions from the r- and s-processes, is moderately 
deficient, and Eu, which has a minimal s-process contri- 
bution, is not deficient at all. The fact that the heavy ele- 
ments in both 16 Cyg A and B show the same abundance 
pattern as a function of Tq suggests that this is a pri- 
mordial effect — the result of Galactic chemical evolution 
(GCE). The slightly super-solar metallicity of the system 
is not in contradiction, as the slope of [ Ba/H] versus time 
is ste eper than that for [Fe/H] (e.g., lEdvardsson et al.l 
[T99I . 

The works by M09 and R09 concentrated only on so- 
lar twin stars, therefore minimizing these relatively small 
GCE effects. Nevertheless, it is interesting to note that 
the solar twin abundances of Y, Zr, and Ba by R09 show 
larger star-to-star error than those of lighter elements 
(see their Figs. 1 and 3), a result that could be associ- 
ated to the age span of their sample. The opposite ef- 
fect of that observed in the 16 Cyg system (i.e., neutron 
capture element abundances higher than the mean solar 
twin trend) has been recently observed in the younger 
solar twin star 18Sco (Melendez et al., in preparation), 
reinforcing the idea that this is related to GCE. Thus, 
our results for the heavy elements tell us that extreme 
caution must be practiced when interpretations of abun- 
dance ratios versus condensation temperature are made 
for stars that are not so similar to the Sun, in particu- 
lar stars that are significantly older or younger, because 
GCE effects are not entirely negligible. 

5.3. Did 16 Cyg A accrete planetary material? Clues 
from lithium abundances 

Based on the discussion given in the previous sections, 
we are led to conclude that the metallicity difference be- 
tween 16 Cyg A and B is not primordial but acquired. 
Either the photosphere of 16 Cyg A became more metal- 
rich or that of 16 Cyg B became more metal-poor after 
they were formed (or while they were being formed) from 
the same molec ular cloud with the sarn e initial chemi- 
cal composition. iLaws fc Gonzalez! (|2001l ) suggested that 
16 Cyg A accreted either 2.5 Earth masses of chondritic 
material or 0.3 Mj of gas giant material to explain its 
higher [Fe/H] (they derived A[Fe/H] = -)-0.025 ± 0.009). 
They also argued that the late accretion of planetary 
material could increase the lithium abundance on the 
star's convective envelope. Indeed, inspection of the 
6708 A region clearly shows that 16 Cyg B is lithium 
poor compared to 16 Cyg A (Fig. jTOj). Our derived Li 
abundances are: ^(Li) = 1.34 ± 0.04 for 16 Cyg A and 
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Figure 10. Spectral region around the Li doublet feature 
(6708 A). Our solar spectrum (solid line) is shown along spectra 
of 16CygA (dotted line) and 16CygB (dashed line). 

^(Li) = 0.73 ± 0.06 for 16CygB, in good agreement 
with previously published estimates ("King et al.' 119971: 
Luck fc Hciter 2006; Taked a et al. 2007; Israelia n et al.l 
20091 : iBaumann et al.l[20l0t ). Thus, 16CygA has about 
4.1 times (~ 0.6 dex) more lithium than 16CygB in its 
photosphere. 

The behavior of Li abundances in solar-type stars 
is very complex and not yet full y understood (e.g., 
iChaboveil [19981 lPinsonneau'rai2010t ). In general, how- 
ever, inspection of large samples clearly shows that at 
around I Mq and Tes — 5780 K a very steep slope of 
increasing Li abundances with stellar mass (or Tes) i s 
present (e.g.. iTakeda et al.|[2007l: IBaumann et al.ll2Q10D . 
Fitting a straight line to the Li abundance versus Tog 
relation of stars in the 5700-5900 K range we find that 
the lithium abundances of stars with Tcff values similar 
to those of 16CygA and B should differ by a factor of 
3.5 (~ 0.55 dex), i.e., very close to the observed value. 
Fish et al. (in preparation) have recently derived lithium 
abundances of more than 700 solar-type stars and com- 
piled and homogenized lithium abundance data for more 
than 1400 FGK stars. They also derive masses and ages 
in a consistent manner employing the exact same method 
used in this work. The lithium versus age trends for 
stars of mass and metal licity simila r to those of 16 Cyg A 
and B, taken from the ' Fish et al.l catalog, are shown in 
Fig. [TTJ The lithium abundances of both 16 Cyg A and 
B are perfectly normal for stars of their mass, age, and 
metallicity, despite the somewhat large scatter, in par- 
ticular for stars like 16 Cyg B. 

If the convective envelope of 16 Cyg A, after it reached 
its present size, accreted 0.3 Afj of planetary material 
tha t has not experien c ed lith ium depletion, as suggested 
by iLaws fc Gonzaled (|2001[ ). the photospheric lithium 
abundance would have increased by about 0.4 dex (a fac- 
tor of 2.5). Thus, if 16 Cyg A was polluted by planetary 
material, considering the mass effect on lithium deple- 
tion, the lithium abundance of 16 Cyg A should be about 
one order of magnitude larger than that of 16 Cyg B, 
which is significantly larger than observed. Such large 
lithium abundance for 16 Cyg A would also make the star 
somewhat lithium rich compared to stars of similar mass 
and metallicity (top panel in Fig. lTlT) . Already the mass 
effect explains most, if not all, of the observed difference. 
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Figure 11. Evolution of lithium abundance for stars of mass and 
metallicity similar to those of 16 Cyg A (top panel) and 16 Cyg B 
(bottom panel). In each panel the stellar mass is within 0.03 M0 
and metallicity within 0.05 dex from those of the 16 Cyg stars. 
Only stars with age/error > 3 are plotted. Upper limits on Li 
abundances are shown with downwards arrows. Open (filled) cir- 
cles represent stars without (with) gas giant planets detected while 
open squares correspond to stars for which planet information is 
not available. Five-pointed stars show the location of 16 Cyg A 
(top panel) and 16 Cyg B (bottom panel). 



The work bv lDelivannis et ahl (J200CO on beryllium abun- 
dances also provides no evidence for accretion of plan- 
etary material by 16 Cyg A, although the authors warm 
that the errors are too large to make a firm conclusion. 

5.4. A signature of gas giant planet formation? 

A different possibility to explain the elemental abun- 
dance differences in the 16 Cyg system is that it was pro- 
duced by the formation of the planet around 16 Cyg B. 
Similar to M09, we propose that the convective zone of 
16 Cyg B has been depleted of metals due to the forma- 
tion of a planet, i.e., that the missing mass of metals in 
16 Cyg B is currently locked-up in its planet. 

Before exploring this hypothesis in detail, we should 
re-address the problem of lithium abundance. A num- 
ber of authors have suggested that lower lithium abun- 
dances could be correlated with the presence of gas gi- 
ants (e.g., Israelian et al. 2009; Gonzalez et al. 2010a), 
although their preferred explanation is that it relates 
to physical processes that affect the rotational evolu- 
tion of stars with planets and hence their lithium de- 
pletion history. A few planet-hosts and single stars are 
plotted in Fig. [TTJ Age, mass, and metallicity effects 
fully explain the observed abundances of both types of 
stars, without invoking the presence or absence of gas 
giant planets, consiste nt wi th the works by, for example, 
IBaumann et al.l ()2010D and lGhezzi et al.l pOlOO . This is 
the case also for the 16 Cyg stars; i.e., the significantly 
lower lithium abundance of 16 Cyg B (~ 0.6 dex relative 
to 16 Cyg A) is not due to the presence of its planet. In 
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our hypothesis, however, a very small fraction of this dif- 
ference (~ 0.04 dex) could be due to the formation of the 
gas giant planet around 16CygB. However, the effect is 
not related to enhanced lithium depletion in the photo- 
sphere of 16 CygB but to the fact that the planet around 
16 Cyg B retained a certain amount of lithium that would 
have otherwise been accreted by the star. 

As shown in Fig. |9l the 16CygA-B abundance dif- 
ferences do not correlate with dust condensation tem- 
perature as in M09. If we force the volatiles to match 
the mean abundance pattern of solar twins, the refrac- 
tories are too low by 0.06 dex and vice- versa (see solid 
lines in Fig. ^. If the metallicity difference between 
16 Cyg A and B is due to the formation of the planet 
around 16 CygB, we must assume that the abundance 
ratios of metals in this planet are very similar to those in 
the star 16 CygB. In particular, we have to assume that 
the volatile to refractory abundance ratios are equal to 
those in the photosphere of 16 CygB (excluding H and 
He, of course). Thus, the planet must be a gas giant, as 
indeed observed. 

We explore our hypothesis quantitatively with a toy 
model. If the planet around 16 Cyg B were to be added 
to the convective zone of the star, the stellar metallicity 
would change by: 



A[M/H] = log 



{Z/X)czMcz + {Z/X)pMp 



(Z/X)czAf, 



cz 



(1) 



where {Z/X)cz is the ratio of the fractional abundance of 
metals (Z) relative to hydrogen (X) in the unperturbed 
convective zone (loosely called "metallicity" in this con- 
text), Mcz the mass of the convective zone, {Z/X)p the 
metallicity of the planet, and Mp the mass of the planet. 
According to lAsplund eFall (|2059l ). (Z/X)'^^ = 0.018. 
Thus, for 16 CygB ([M/H]=0.061) we adopt iZ/X)cz = 
jqO.061 X 0.018 = 0.021. 

We can use Eq. [T] to calculate the mass of the convec- 
tive envelope necessary to explain the observed metal- 
licity difference of A[M/H] = -1-0.04 given the plane- 
tary parameters {Z/X)p and Mp. Since a minimum 
mass estimate of l.S Mj is available from the work by 
iCochran et al.l (|1997t l , we explore Eq. [1] for a range of 
planet masses from 1.5 to 9.5 M j. Th e results from our 
toy model are illustrated in Fig. I12l^'^l 

This toy model is intended to counteract the immedi- 
ate aftermath of the planet formation process. Thus, the 
mass of the convective envelope that appears in Eq. [1] 
corresponds to the convective zone at the time of planet 
formation. Standard models of early stellar evolution 
suggest that solar analog stars begin fully convective and 
only after about 30 Myr their co nvective envelopes reach 
their present-d ay masses (e.g., ID'Antona fc Mazzitellil 
119941: iSerenelhet al...20ll[ ). On the other hand, obser- 
vations of disks around young stars suggest a shorter 
timescale for the process of planet formation, of order 
10 Myr (e.g., iMamaiekl [20091: iMeveii [2009t iFedele et al.l 

^* Our simple model assumes instantaneous accretion of metal- 
deficient material, therefore implying that the planet itself formed 
instantaneously, which is a very crude approximation. More de- 
tailed models should consider the finite formation time scale of the 
planet in addition to the temporal evolution of the mass of the 
convective envelope but that is beyond the scope of our work. Our 
results are therefore only useful as first order approximations. 
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Figure 12. Solid lines correspond to the mass of the convective 
envelope of 16 Cyg B necessary to explain the metallicity difference 
between 16 Cyg A and B as a function of the metallicity of the 
planet around 16 CygB and for planet masses between 1.5 and 
9.5 Mj in increments of 1.0 Mj. Dashed lines correspond to the 
mass of the convective envelope of 16 Cyg B for stellar ages between 
5 and > 30 Myr according to the standard model bv lSerenelli et al.l 
II2011I ). Dotted lines represent a similar evolution for Mqi^ from 
to 10 Myr, but according to one of t he non-stan dard models with 
episodic accretion by Baraffe & Cha brierl 1)20101 ). 



I2010f ). Thus, we investigate our simple model keeping 
Mcz as a variable quantity but const raining it s tem- 
poral evolution as computed by Seren elli et al.l (|201lL 
their Fig. 1). The latter corresponds to a IMq star of 
solar-metallicity. To take into account the small differ- 
ences in mass and [ Fe/H] for 16 CygB, we scaled the 
iSerenelli et al.l (|2011l ) result using the present-day con- 
vective envelope m a sses a s a function of Toff, [Fe/H] by 
iPinsonnea ult et aD (|2001[ ). The dashed lines plotted in 
Fig. [12] correspond to the mass of the convective zone at 
different stellar ages according to standard models. In 
this figure we also show, with dotted lines, the time evo- 
lution of Mqz for o ne of the non-standard mo dels with 
episodic accretion bv lBaraffe fc Chabrieii ()2010D . We dis- 
cuss this case later in this section. 

According to our pollution model with standard Mqi 
evolution, if the planet has (Z/X)p = 0.1 and Mp = 
1.5 Mj, it should have formed when the star was about 
25 Myr old but if the planet has Mp = 9.5Mj, it should 
have formed 10 Myr after the star was born, when the 
convective envelope was more massive, therefore dilut- 
ing more the larger amount of metal-depleted gas ac- 
creted. The planet around 16 CygB could not have 
formed 25 Myr after the stellar birth since the convec- 
tive envelope at that point was too small. Only if the 
metallicity of the planet is significantly lower (< 0.03) it 
could have formed at these late times. On the other 
hand, if the planet is more metal rich, for example 
{Z/X)p = 0.7, it must have formed in the first 10 Myr 
since the birth of the star, and earlier for larger planet 
masses. Such large planet metallicity appears, how- 
ever, unphysical for gas giants. The amount of met- 
als in Jupiter (cor e and envelope) is abo ut 12-37 Mq 
(jGuillot 2005: Fortnev fc Nettelmannll2010f ) and it seems 
that all giant planets have at least 10-15 M^ of met- 
als (jMiller fc FortnevI 1201 ID . Since Jupiter's mass is 
318 M0, a total amount of metals of 12-37 M^ would 
imply {Z/X)p = 0.04 - 0.12. 

Theoretical models of planet formation predict that 
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gas giant planets form relativel y quickly, on time 
scales shorter than lOMyr (e.g., iPoUack et al.l 119961 : 
iGuilera et al.ll20l"ll ). Yet, for reasonable values of the 
planet's metallicity and mass, planet formation time 
scales from our pollution model with standard Mcz evo- 
lution are in the 10-30 Myr range, i.e., we require the 
process of gas giant planet formation to be significantly 
longer. The exact same problem is faced by M09 sce- 
nario for the formation of terrestrial planets. They have 
to assume that while the terrestrial planets were being 
formed, the mass of the solar convective envelope had 
already reached its present size. A possible solution to 
this problem may lie in the new developments in the field 
of early stellar evolution. In particular, the impact of 
episodic accretion on the internal structure of stars may 
help solving this discrepancy. 

Models for the early evolution of stars with non-steady 
accretion have been computed by Baraffe k, Chabrierl 
()2010[) . They demonstrate that the internal evolution 
of stars is highly dependent on the accretion history. Of 
particular interest for our work is the fact that their ac- 
cretion models predict higher central temperatures and 
hotter overall structures leading to a more rapid forma- 
tion of the radiative core in solar-type stars. In one of 
their models, the mass of the convective envelope of a 
1 Mq star reaches its present-day size after only about 
10 Myr from the stellar birthline. The evolution of Afcz 
for this model in particular is represented with the dotted 
lines in Fig. [T2l 

One possibility to explain the metallicity difference be- 
tween 16Cyg A and B in the context of gas giant planet 
formation, and at the same time remain compatible with 
the terrestrial planet signature hypothesis by M09, con- 
sists in assuming that the gas giants formed during the 
first ~ 5 Myr after the birth of the star whereas the ter- 
restrial planets formed later (> 10 Myr). More specifi- 
cally, the timescales of accretion of metal-deficient and 
fractionated gas should be consistent with this assump- 
tion. In addition, the evolution of Mcz should be quicker 
than predicted by standard models; i.e., Afcz should 
reach its present d ay mass in about 10 Myr, a s it is the 
case of some of the iBaraffe fc Chabried (|2010( ) models. 

5.5. Elemental abundances in the solar nebula 

An interesting consequence of the proposed scenario 
is that the solar nebula would have been more metal- 
rich than the present-day solar photosphere. In fact, we 
could re-construct the original solar abundances if we had 
a good knowledge of the composition of the gas giants. 
The combined mass of the gas giants in the solar system 
is about 1.4 Af j, slightly lower than the minimum mass of 
the planet around 16CygB. On the other hand, the con- 
vective envelope of the Sun is somewhat less massive than 
that of 16 Cyg B. Given the uncertainties in the many pa- 
rameters involved to calculate the metallicity shift due to 
the formation of the gas giants in the solar system, we 
will assume for simplicity that their impact was simi- 
lar to the case of 16 Cyg B. Thus, the formation of the 
gas giants in the solar system could have shifted the so- 
lar abundances of all elements down by about 0.04 dex. 
Later, the formation of rocky planets lowered even more 
the abundances of refractory elements by up to 0.08 dex 
while keeping the volatile abundances nearly constant 
(cf. M09, R09). Thus, for example, the abundances 



of some important elements in the solar nebula would 
have been [C,N,0/H] ~ -1-0.05, [Fe/H] ~ -1-0.09, and 
[Al, Sc,Ti/H] ~ 0.12. Note that in this calculation we 
have ignored the impact of diffusion, which would fur- 
ther decrease the present-da y observed sola r abundances 
by ab out 0.05 dex (Turcotte &: Wimmer-Schweingrubeil 
|2002|) . Naturally, if the mass of the planet around 
16 Cyg B is significantly larger than 1.5 Af j, the impact of 
the formation of Jupiter on the solar abundances would 
be lower than 0.04 dex. 

The elemental abundances measured in CM chondrite 
meteorites, which are often assumed to have retained the 
solar nebula composition, agree very well with th ose mea- 
sured in the present-day so lar photosphere ('e.g.. lLoddersl 
120031 lAsplundet aIll2009D . The impact of planet for- 
mation and diffusion on the solar abundances described 
above is not necessarily in contradiction with this ob- 
servation because when these comparisons are made the 
meteoritic abundances have to be normalized using the 
solar abundance of Si (or a mixture of several metals 
with low abundance uncertainties), therefore removing 
any overall metallicity offsets. The trend with condensa- 
tion temperature, however, would not disappear in this 
way but it has not yet been detected, probably owing 
to the still relatively large errors in the determination of 
absolute solar elemental abundances. 

According to our proposed scenario, the differences be- 
tween the initial and present-day solar abundances are 
not negligible and they could be relevant for example in 
studies of Galactic chemical evolution, chemical tagging 
of dynamical groups, and stellar interior modeling. Fur- 
ther work on the impact of planet formation on stellar 
abundances, both from the observational and theoretical 
sides, is therefore urgently needed. 

6. SUMMARY 

We have determined the stellar parameters and abun- 
dances of 25 elements in 16 Cyg A and 16 Cyg B with very 
high precision. Their lithium abundances appear per- 
fectly normal considering the observed lithium depletion 
patterns in field solar- type stars as a function of mass and 
metallicity and the small (but non-negligible) differences 
in these two parameters for the 16 Cyg stars. Accretion 
of planetary material onto 16 Cyg A and enhanced pho- 
tospheric lithium depletion due to the presence of the 
planet around 16 Cyg B are not necessary to explain the 
observed lithium abundances. 

Excluding lithium, we show that 16 Cyg A is more 
metal-rich than 16 Cyg B by A[M/H] = -f 0.041 ± 
0.007 dex. The abundance differences detected do not 
correlate with dust condensation temperature. Thus, if 
the difference is due to the accretion of metal-depleted 
gas onto 16 Cyg B, the abundance ratios of different met- 
als in this gas, in particular the volatile to refractory 
abundance ratio, should be nearly the same as that found 
today in the photosphere of 16 CygB. We propose a sce- 
nario in which both 16 Cyg A and B started with the 
exact same composition but the formation of a planet 
around 16 CygB made this secondary star slightly metal- 
poor. The composition of the planet must be the same 
as that of the stars (with the exception of H and He) in 
order for this hypothesis to work on a quantitative basis. 
Moreover, the planet must have formed while the con- 
vective envelope of the star was still relatively massive 
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(about 10% of the stellar mass). 

The idea proposed to explain the abundance differences 
in the 16 Cyg binary syst em is similar to that suggested 
by iMelendez et~all (|2009[ ) and iRamirez et all (|2009t ) to 
explain the anomalous solar abundances as a signature 
of the formation of terrestrial planets. The main differ- 
ence is that the planet around 16CygB is a gas giant. 
We can reconcile these two hypothesis by assuming that 
gas giants form while the convective envelope is massive 
whereas terrestrial planets form later, when the convec- 
tive zone reaches its present-day size. In any case, the 
formation of the stellar radiative core should be quicker 
than that predicted by standard stellar evolution mod- 
els. The impact of episodic accretion on early stellar 
evolution could provide a solution to this problem. 

If this scenario is correct, we must conclude that the 
Sun's original metallicity was higher, an effect possibly 
enhanced by atomic diffusion. Perhaps the regions im- 
mediately below the base of the solar convective enve- 
lope are significantly more metal-rich than the photo- 
sphere, a condition that has been proposed as a possi- 
ble solution to the so-called solar abundance crisis (e.g., 
[Nordlund 2009), i.e., the severe discrepancy between the 
observed photospheric solar abundances and those in- 
ferred from h elioseismological dat a and s tellar interior 
model s Ce.g.. iBasu fc Antial 12004 l2008t iBahcall et all 
120051: iDelahave fc Pinsonneaulti 120061 ). Note, however, 
that this idea seems not to work very well on a detailed, 
quantitative basis (e.g.. lSerenelli et al.ll201ll ). 

Our detail study of the 16 Cyg system shows that high 
precision elemental abundance work in wide binary stars 
with planets can provide us with important clues to un- 
derstand the impact of planet formation on the chemical 
composition of stars. We encourage future spectroscopic 
studies of binary stars with similar components, both to 
measure small differences in their chemical composition 
and to look for planets, especially if the stars are solar 
analogs. 
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